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Abstract 

This paper discusses the dimension of spHne spaces S(m, n,m — l,n — l, 
over certain type of hierarchical T-meshes. The major step is to set up a bijec- 
tion between the sphne space S(m, n, m — 1, n — 1, and a univariate sphne 
space whose definition depends on the 1-edges of the extended T-mesh. We 
decompose the univariate spline space into direct sums in the sense of isomor- 
phism using the theory of the short exact sequence in homological algebra. 
According to the decomposition of the univariate spline space, the dimension 
formula of the spline space S{m,n,m — l,n — 1, J7) over certain type of hi- 
erarchical T-mesh is presented. A set of basis functions of the spline space is 
also constructed. 

Keywords: Dimension formula, spline space, T-mesh, homology. 



1 Introduction 

Non-Uniform Rational B-Splines (NURBS) are popular tool to represent surface 
models in Computer Aided Geometric Design (CAGD) and Computer Graphics. 
However, due to the tensor-product structure of NURBS, local refinement of surface 
models based on NURBS is impossible, and NURBS models generally contain large 
number of superfluous control points. To overcome the above drawbacks, Sederberg 
et al. introduced T-splines, the control meshes of which allow T-junctions ([H [2])- 
T-splines provide local refinement strategy and can reduce the large number of 
superfluous control points in NURBS models. 
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In [3], the concept of spline spaces over T- meshes is introduced. Different from T- 
sphnes, a sphne over a T-mesh is a single polynomial within each cell of the T-mesh, 
and it achieves the specified smoothness across the common edges. Spline spaces 
over T-meshes are suitable for geometric modelingj4| and they can be applied for 
analysis naturally, since it is easy to do standard Finite Element Method analysis 
based on splines over T-meshes [5]. 

One major issue in the theory of splines over T-meshes is to study the dimension 
of the spline spaces. There have been a few literature focusing on the problem so far. 
In [3], the dimension formula for the spline space S{m,n,a, (3, ^) is obtained with 
constraints m > 2a+l, n > 2/3+1. The result is further improved by Li et al. [6]. For 
the spline spaces with highest order smoothness where m < 2a, n < 2/3, the authors 
of the current paper derived a dimension formula for biquadratic spline spaces 
(that is, S(2,2, 1, 1, ^)) over hierarchical T-meshes[7]. Recently, B. Mourrain gives 
a general formula for the spline spaces S(m, n, a, (3, ^) by homological techniques [5]. 
Unfortunately there is a term in the dimension formula which is very hard to compute 
in practice. 

Contrary to the above positive results, Li and Chen have showed that the di- 
mension of spline spaces with highest order smoothness over T-meshes may depend 
on the geometry besides the topology information of the T-meshes [9]. The result 
suggests that it is vain to study the dimension formula over general T-meshes in the 
case of highest order of smoothness. In this paper, the dimension formula of spline 
spaces S(m, n, m — 1,?t, — 1, ^) over certain type of hierarchical T-meshes will be 
explored. 

There are several methods for establishing the dimension of spline spaces, such 
as the B-net method [3], and the smoothing cofactor method [10] and homology 
method [11]. In the smoothing cofactor method, the cofactor of a spline associated 
with the common edge of two adjacent cells is a univariate polynomial. This is 
similar to homology theory in topology. Therefore, using the smoothing cofactor 
method and according to a result in [7], we construct an isomorphism between a 
spline space with highest order smoothness over a T-mesh and a univariate spline 
space satisfying some conditions. For a certain type of hierarchical T-mesh, by giving 
an order of the interior l-edges of the extended T-mesh, we can then decompose the 
univariate spline space into a direct sum, and thus give a dimension formula for the 
spline space over hierarchial T-meshes. 

The rest of the paper is organized as follows. In Section 2, the definitions and 
some results regarding T-meshes and spline spaces over T-meshes are reviewed. In 
Section 3, an equivalence is set up between the spline over a T-mesh and a univariate 
spline space. The proof of the dimension formula is presented in Section 4. Some 
examples are also provided. In Section 5, we conclude the paper with future research 
problems. We leave the proof for a key lemma and the construction of a set of basis 
functions of the spline space in the appendix. 
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2 Spline Spaces over T-meshes 



In this section, we first review some concepts about T-meshes and sphne spaces over 
T-meshes. 



2.1 Spline spaces over T-meshes 



A T-mesh is a rectangular grid that allows T-junctions. For simplicity, in this paper 
we consider only regular T-meshes whose boundary grid lines form a rectangle. We 
adopt the same definitions for vertex, edge, and cell as in |3], and the definitions 
for 1-edge, interior 1-edge, associated tensor product mesh are borrowed 
from [7j. A grid point in a T-mesh is called a vertex of the T-mesh. Vertices 
of a T-mesh are divided into different types. For example, in Figure [II {6j}^Z}° 
are boundary vertices and {fjj^^f are interior vertices. V2 is a crossing vertex and 
{^i}i=i'' U {t^j}i=i — {"^2, &3, ^s} are T-vertices. The line segment connecting two 
adjacent vertices on a grid line is called an edge of T-mesh such as f4f5, &9&10, "^^2^3 
in Figure [TJ h2V2, is a large edge (1-edge for short), which is the longest possible 
line segment consisting of several edges. The boundary of a regular T-mesh consists 
of four l-edges,which are called boundary l-edges. The other l-edges in T-mesh are 
called interior l-edges. A regular T-mesh can be extended to a tensor product mesh, 
called the associated tensor-product mesh, by extending all the interior large edges 
to the boundary, see Figure [H 
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Figure 1: A T-mesh 3^ and its associated tensor product mesh S^'^. 



Given a T-mesh ^ is the set of all the cells of and is the region occupied 
by cells in ^ . Spline spaces over T-meshes are defined by 

S(m, n, a, /3, ^) : = {/(a;, y) e C"'^(n) ■.f{x,y)\4>e P„„, V0 G ^}, (1) 

where P^^ is the space of all the polynomials with bi-degree (m, n), and C"'^ is the 
space consisting of all the bivariate functions that are continuous in Q with order a 
along the x direction and order /3 along the y direction. In this paper, we will focus 
on the spline space S(m, n,m — l,n — 1, 
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2.2 Hierarchical T-meshes and extended T-meshes 

A hierarchical T-mesh is a special type of T-mesh that has a natural level struc- 
ture [1] . It is defined in a recursive fashion. Initially a tensor product mesh (level 0) 
is presumed. From level k to k + 1, a. cell is subdivided at level k into four sub-cells, 
which are cells at level k + 1, by connecting the middle points of the opposite edges 
with two straight lines. Figure |2] illustrates a sequence of hierarchical T-meshes. 



Level Level 1 Level 2 

Figure 2: Hierarchical T-meshes. 

For a T-mesh 3^, the extended T-mesh associated with ^ is an enlarged 
T-mesh by copying each horizontal boundary line of ^ m times, and each vertical 
boundary line of n times, and by extending all the line segments with an end 
point on the boundary of ^ |7]. This can be made precise by the following example. 
Figure [3] illustrates a T-mesh (left) and the extended T-mesh (right) associated with 
degree (3, 3). 



Figure 3: A T-mesh 3^ and its extended T-mesh associated with degree (3,3). 
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2.3 Homogeneous boundary conditions 

A spline space over a given T-mesh with homogeneous boundary conditions is 
defined by 

S(m, n, a, /3, ^) := {/(x, G C"'^(M2) ; /(a;, |^ g P^„, V0 G ^ and = 0}, 

where P^n,^ is defined as before. One important observation in [7] is that the 
two sphne spaces S(m, n, m — 1, n — 1, ^) and S(m, n, m — 1, n — 1, ^'^) are closely 
related. 

Theorem 2.1. |^ Given a T-mesh , and let be the extended T-mesh associated 
with S{m, 71,171 — l,n — 1, . Then 

S(m, n,m — l,n — 1, = S(m, n,m — l,n — 1, (2) 
dim S(m, n, m — 1, n — 1, ^) = dim S(m, n, m — 1, n — 1, (3) 

Based on the above theorem, we only have to consider spline spaces over T-meshes 
with homogeneous boundary conditions. 

3 Equivalent spline spaces 

In this section, we will give an equivalent description of spline spaces S{m,n,m — 
1, n — 1, ^'^) by using the smoothing cofactor method. 

3.1 Smoothing cofactor method 

As before, let ^ and be a T-mesh and its extension respectively. Let ^ = 
{Ci, C2, ■ ■ ■ , Ck} be the set of all the cells in and fli be the region occupied by 
d G ^. Denote fifc+i = U^L^ and = {lli, (^2, ■ ■ ■ , ^^fc+i}- 

Let Ui, U2, U3, G be four regions in adjacent positions as shown in FigureHl 
For a spline function f{x,y) G S(m, n, m — l,n — 1, denote fi{x,y) be the 
bivariate polynomial which coincides with f{x,y) on Ui, i = 1,2,3,4. Note that if 
fi{x,y) = fj{x,y) for two adjacent cells Ui and Uj, then Ui and Uj can be merged 
into a single region, and in this case (xq, yo) is a T-junction. 

By the smoothing cofactor method, we have the following relationship between 
Mx,y): 

Lemma 3.1. JE, \^ Let fi{x, y),i = 1,2,3,4 be defined as in the preceding paragraph. 
Then there exist a constant A; G M and polynomials a{y) G Pn[|/]; b{x) G Pm[a;] such 
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Figure 4: Smoothing conditions in adjacent cells 



that 

fiix,y) = f2ix,y) + h{x){y-yQY, 
h{.x,y) = f2ix,y) + a{y){x - Xq)"", 

U{x, y) = f2{x,y) + a{y){x - xq)™ + b{x){y - yoT + - ^^Thj - VoT, 

where {xo,yo) is the vertex of 3^"^ in Figure\^ Furthermore, a{y),b{x) and k are 
uniquely determined by fi{x,y), i = 1,2,3,4, specifically 

. f d'^h{xo,y) _ d"'f2{xo,y) \ 

""^y' m\ V dx-^ ) ' 

^^""^ ~ n\\ dy- dy- ) ' 

and 

1 l d^+^{h{x,y) + U{x,y)-h{x,y)-h{x,y)) 
ml n\ dx"^dy'^ 

a{y) and b{x) are smoothing co-factors associated the edges between adjacent cells. 
The constant k is called the conformality factor associated with the common vertex 

Based on the above lemma, we convert the dimension problem into the study of 
a univariate spline space. We first introduce the following definition. 

Definition 3.1. Let be a horizontal l-edge in S/'^ , and {f i, f2, ■ ■ ■ , f^} be vertices 
on E^. Assume the x-coordinates of the vertices are Xi < X2 < ■ ■ ■ < Xr respectively. 
We define a univariate spline space with homogeneous boundary conditions associ- 
ated with E^: 

S(m,m - 1, E"^) := {p{x) G C'^-\R) : p{x%,,^^^,^ = Pi{x) G P„[x], 

i = 1,2, . . . ,r - 1, one? p(x)|iR\[^.^,^^] = 0} . (7) 
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Similarly, for a vertical l-edge E"" , we can define a univariate spline space S(n, n — 
1,E^) associated with it. 

By the homogeneous boundary conditions, it is easy to see that for any p{x) G 
S(m, m — 1, E^), there exist constants ki, . . . ,kr such that 



i=l 



Equation is equivalent to a hnear system associated with E^: 

e:=i h = 0, 

kiXi 0, 

' ' ' ) 

I e:=i = 0. 

For a vertical l-edge E^, there is a similar equation 

s 

J2k,{y-y,r^0, 

i=l 

which is also equivalent to a linear system associated with E"": 



(9) 



(10) 



r e:=i h = 0, 

Ei=l ^iVi = 0' 
' ' ' ) 

I ELi hy? = 0. 

By iQ and ffTTj) . one immediately has 
Lemma 3.2. 

dim S(m, m — 1, .E^) = {r — m — 1)+, dim S(n, n — 1, 
ifere ti+ = max(0,'u). 



(s - n- 1)+. 



3.2 Conformality vector spaces 

To study the dimension of spline space S(m, n,m — l,n — l, ^^), we need to consider 
the conformality conditions for all the (horizontal and vertical) l-edges. Thus we 
introduce the following definition. 
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Definition 3.2. Let , i = 1,2, ... ,p he all the horizontal l-edges of , and E^, 
j = 1,2, ... ,q be all the vertical l-edges of S/''^ . Define a linear space W[3/''^] by 



W[^'] := {k={h,h 



■.L 



0, i 



,P, 3 



where v is the number vertices of fcj is the conformality factor corresponding to 
i-th vertex of 3"^ , and = and = are linear systems associated with the 
l-edges E^ and Ej respectively. W[3/''^] is called the conformality vector space of 
S{m,n,m — l,n — 1, ^^). Similarly, one can define the conformality vector space 
W[E] of S{m,m — 1, E^) (or S{n,n — 1, E^)) associated with a l-edge. 

The following facts should be noted regarding W[^'^]. 

1. W\S/'\ can also be defined over a general T-mesh ^ (besides an extended 
T-mesh). For any f{x, y) G S(m, n, m — 1, n — 1, there is a unique corre- 
sponding vector k G W[^] called conformality vector. 

2. ki, which is the i—th component of k, is the conformality factor corresponding 
to i—th vertex in JT'^. This vertex is the intersection point of two l-edges. 
So ki has to satisfy both equations Qi and (fTTl) associated with the two 1- 
edges. Once k is determined, the smoothing co-factors a{y) and b{x) can 
be constructed accordingly. Figure [5] illustrates an example, where a{y) = 

h{y - yiT + ^4(2/ - 2/2)" + h{y - vsT, K^) = h{x - xi)™ + k2{x - X2T'. 

a{y) is determined by k^, k^, k^ whose corresponding vertices lie beneath the 
edge associating with a{y); b{x) is determined by ki,k2 whose corresponding 
vertices lie on the left of the edge corresponding to b{x). 
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Figure 5: Conformality vector 



3. A naught value of ki which is the conformality factor corresponding to a T- 
vertex results in the vanishing of the corresponding vertex in the T-mesh 
£/' . An example is illustrated in Figure |6] where k-j = 0, and in this case 
S(m, n,m — l,n — 1, ^) = S(m, n,m — l,n — 1, ^'), or ^ degenerates to 
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Figure 6: The naught value of a conformahty factor. 



As an example, we discuss the conformality vectors of B-spline functions. 

Example 3.1. Let A : — oo < Xi < X2 < ■■■ < Xr < 00 (r > m + 1) be a 
partition ofW, and S(m, m — 1, A) be the B-spline function space defined over A, 
that is, S{m,m — 1,A) consists of piecewise polynomials of degree m with C"^~^ 
continuity over A. The B-spine basis function iV[a;i, Xj+i, . . . , Xj+m+i] is an ele- 
ment in S{m,m — l,Ei), where Ei is the interval [xj, Xj+m+i] with the partition: 
Xi < Xj+i < . . . < Xj+m+i ■ Here l<i<r — m — 1. The conformality vector 
k = {ki, fcj+i, . . . , fcj+m+i) of N[xi, Xi+i, . . . , Xi+m+i] Can be obtained by solving the 
associated linear system 1^. It is easy to see that kjkj^i <0, j = i,i-\-l,...,i-\-m 
and kj > 0. 



Example 3.2. Let S(m, m — 1, A^;) be the B-spline function space defined over 
Ax : —00 < Xi < X2 < ... < Xr < 00 (r > m -\- 1), and S{n, n — 1, Ay) be 
the B-spline function space defined over Ay : —00 < j/i < ?/2 < • • • < 2/s < 00 
(s > n -\- 1). We denote as a tensor product mesh A^ x Ay. It is easy to see 
that, /(x) G S(m, m — 1, A^.) and g{y) G S{n,n — I, Ay) implies that f{x)g{y) G 
S(m, n, m — l,n — 

Now assume that the B-spline basis functions N[xi, Xj+i, . . . , Xi+m+i] G S(m, m — 
1, Ax) and N[yj,yj^i, . . . , yj^n+i] ^ S{n, n — 1, Ay) have conformality vectors : = 



j+i' 



1 '^i+m+lj 



and k 



^2 



respectively. Then by equa- 

, -^i+m+l] ■ A^bj! 2/i+l! • • • ; Vj+n+l] ^ 



{kl,k, 

tion (Q, the B-spline basis function N[xi, Xj+i, . . . , x 
S{m,n,m — l,n — 1, ^) has a conformality vector ]<} ® k^ which is a vector of 
dimension {m + 2)(n + 2) with elements k^k"^, p = i,i -\- 1, . . . ,i -\- m -\- 1, q = 
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3.3 Equivalence of spline spaces 

Based on the above preparations, we obtain a mapping between the sphne space 
S(m, n,m — l,n — 1, ,9''^) and the conformahty vector space W[S/''^\. 

jr : S(m,n,m-l,n-l,e^^) — >W[3r']. (12) 

By Equation ([6]), is a hnear mapping due to the hnear property of the operator 
Qm+n jQx^dif. In fact, is an isomorphic mapping. 

Theorem 3.3. The mapping : S(m, n,m — l,n — l, — > W[3^^] is bijective. 

Proof. We first prove that is injective. It is enough to show that Jt^f = implies 
/ = for / G S(m, n, m — 1, n — 1, ^^). By the remark 3 following Definition 13. 21 if 
the conformahty vector of f{x, y) is a zero vector, then f{x, y) is a single polynomial 
over By the homogeneous boundary conditions, f{x,y) = 0. Thus J(f is 

injective. 

Next we show is surjective. By the remark 2 following Definition 13.21 for a 
given conformahty vector k G one can construct a smoothing cofactor for 

each edge of the T-mesh and thus obtains a spline function f{x, y) G S(m, n, m— 
l,n — 1, ^'^) corresponding to k, that is, J(f is surjective. Thus the mapping J(f is 
bijective. □ 

By the above theorem, the spline space S{m,n,m — l,n — 1, is isomorphic 
to the conformahty vector space 11^ [^^], and 

dimS(m,n,m- l,n- 1,^^) = dimW^[^"]. (13) 

Similarly, one can show that 

S{m,m-1,E'') = W[E''], S{n,n - 1, E'') = (14) 

for a horizontal 1-edge and a vertical l-edge E"". 

In the following, we only have to analyze the structure of W[^^]. 

4 The dimension formula 

In this section, we will derive a dimension formula for the spline S(m, n,m — l,n — 
1, 3^"^) over a certain type of hierarchical T-meshes by direct sum decomposition 
of the spline space. The decomposition is according to each interior 1-edge of the 
T-mesh A set of basis functions of the spline space is also constructed. 
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4.1 Some definitions 



We first introduce a certain type of hierarchical T-meshes associated with degree 
m,n and denote it by 3^rn,n- 

Definition 4.1. Let = [xo,xi,--- ,Xp] x [yo,yi,--- ,yq] be a tensor product 
mesh, and let X = {xo,Xm_i,-- - ,Xs{m-i),Xp} and Y = {yo, 1/n-i, • ■ ■ ,yt{n^i),yq}, 
where s and t are unique integers such that < p — s{m — 1) < m — 1,0 < 
q — t{n — 1) <n — l. The domain f2 = [xojXp] x [?/o;Z/g] is subdivided by the lines 
{x = Xi,y = yj : Xi G X,yj G Y} into sub-domains, each of which is occupied by 
a local tensor product mesh of size (m — 1) x — 1) (or smaller near the right 
or upper boundary lines). Each subdomain is called a {m,n)-subdomain of a 
{m,n)-subdomain is called s boundary one, if it is near the boundary lines of T- 
mesh. A subdomain is said subdivided if each cell in the subdomain is subdivided ( a 
cell is subdivided by connecting the middle points of oppose sides of the cell with two 
line segments). A subdomain is called isolated if it is subdivided while its adjacent 
subdomains (two subdomains are called adjacent if they share a common boundary 
line segment) are not subdivided and it is not boundary one. 

Figure [7](a) illustrates an example where a tensor product mesh is subdivided 
into (3, 3)-subdomains by blue lines. In Figure [7](b), two subdomains are further 
subdivided by red lines, and these one subdomain is isolated. 



(a) (b) (c) 

Figure 7: Hierarchical T-mesh ^rn,n in the case of m = n = 3, p = 5, g = 6 



Definition 4.2. A hierarchical T-mesh ^m.n associated with [m, n) is defined level 
by level. At level k = 0, S/'m.n is a tensor product mesh From level to level 1, 
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subdivide ^ into {m,n)-subdomains, and some of which are further subdivided to 
get the mesh at level 1. Generally from level k to level k + 1, subdivide some local 
tensor product meshes at level k (which are obtained by subdividing subdomains at 
level k — 1) into (m, n) -subdomains at level k, some of which are further subdivided 
to get level k + 1 mesh. 

The T-mesh ^rn,n is a special type of hierarchical T-mesh [3]. In particular, it is 
a regular hierarchical T-mesh when m = n = 2. Figure [7] illustrates the process of 
generating a hierarchical T-mesh £^3^3. Figure [Tl^a) shows that a tensor product mesh 
is subdivided into subdomains at level k = 0. Figure [7](b) shows two subdomains 
at level k = are subdivided to get a mesh at level k = 1. Figure [7|(c) shows 
two local tensor product meshes (corresponding to the two subdivided subdomains 
in Figure [71(b)) are subdivided into (3, 3)-subdomains at level k = 1, and three 
subdomains at level k = 1 are further subdivided to get a mesh at level k = 2. 

4.2 The dimension formulas 

Now we are able to state our main result in the paper - the dimension formula 
for spline space S{m,n,m — l,n — 1, ^„), where is =v^,„'s extended T-mesh 
associated with degree m, n. 

Theorem 4.1. Let be the number of crossing vertices of Assume ^„ 

has Eh interior horizontal l-edges and Ey interior vertical l-edges. Then 

dim S(m, n, m — 1, n — 1, ^ ^) 

= (m - l){n - 1) + y+ - (m - 1)Eh - {n - l)Ev + 5, (15) 
where 5 is the number of isolated subdomains of S/^^ „ at all levels. 

Before proving the theorem, we need some preparations. 

Lemma 4.2. Let £^rn,n be a hierarchical T-mesh associated with {m,n), and ^„ be 
its extended T-mesh. E is the set of all the interior l-edges of H^^^. Then the l-edges 
in E can he ordered as Et -< Et-i -<■■■-< Ei such that the projection mapping 

TT : W[^^ ^ W[E,] 

is surjective, z = 1, 2, . . . , t. Here ^ is the T-mesh obtained by deleting {Ei, ■ ■ ■ , Ei_i} 
from S^^n, and Ei d S/'i is regarded as a I- edge in J^. 

Proof. See the appendix. □ 

Base on the above lemma, we can decompose the conformality vector space 
into direct sums of W^f-Ei]. 
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Theorem 4.3. Let the notations he the same as in Lemma 



El 



Then 



w\sr^, ^ w\sr,^r\ © z = 1, 2, . . . , t. 

Proof. Since = ^i^i U S(m, n, m — 1, n — 1, C S(m, n, m — 1, n — 1, .^3^). 

Correspondingly, W[^i+i] can be regarded as a subspace of W[^i\ by taking every 
vector k G W[3^i+i] as a vector in W[^] whose components corresponding to the 
vertices of Ei are zero, and the remaining components are the same as k's. 

Now consider the sequence 

^ -^-^ W[^] W[E,] , (16) 

where i in the natural embedded mapping, and vr is the projection mapping defined 
in Lemma [4.21 We will show that the sequence is exact. Since i is injective and tt is 
surjective by Lemma [4.2[ it is enough to show that Im{i) = Ker{T[). 

On one hand, for any k G W[£^,^i], i(k) = k' is a vector in W[£^] whose com- 
ponents (conformality factors) corresponding to the vertices of Ei are zero. Thus 
7r(i(k)) = G W[Ei], that is, /m(i) C Kerin). 

On the other hand, for any k G Ker{7T), 7r(k) = 0, that is, the components of k 
corresponding to the vertices of Ei are zero. Therefore k G Im{i), i.e., Ker{iT) C 
Im{i). 

Thus the sequence ( IT6ll is exact. Since W[^i+i], W[£^] and are linear 

spaces, the sequence is split, and therefore W[£^] = W[^i+i] © as required. 

□ 

Based on the above theorem, we have 

Corollary 4.4. The set of functions U-^]^{5*(x, ?/)}^~^'~'"~"^ defined in the proof of 
Lemma \4.^ form a basis for the spline space S{m,n,m — l,n — 1, ^„). 



Now we are ready to prove the dimension formula ( JTSll . 



Proof of Theorem \4.1\ By Theorem 14.31 can be decomposed into direct 

sums of 

The interior l-edges E = {Ei,E2, . . . ,Et} of can be divided into disjoint 
sets E\ i = 0,1,..., I, i.e., E = uJ^qE' and E* n E^' = 0, i ^ j. E* is the set of 
interior l-edges of defined at level i, i = 0,1, . . . ,1. In particular, E" is the set 
of the interior l-edges of the initial extended tensor product mesh By the order 
defined in Lemma 14.21 any element in E* precedes any element in E-' for i > j. 
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By the direct sum decomposition of W[^j^^], 

I 

dimW[^lJ = E E d\mW[E,]. (17) 

Now we count dimVr[£'j] according to the level of Ej from the highest level / to 
the lowest level 0. 

By Lemma IX^and ffTI]) . 

dimiy[Ej] = dimS(m,m - 1, Ej) = (f+ - m + 1) + 

if Ej is a horizontal l-edge, and 

dim 11" [Ej] = dim S(ra, n - 1, Ej) = {vf -n + l) + 

if Ej is a vertical l-edge. Here is the number of crossing vertices of Ej in the 
mesh ^j. 

Consider all the interior l-edges E* in a fixed level i {i = 1,1 — 1, . . . ,0). E* is 
obtained by subdividing some subdomains at level i — 1 (except for i = where E° 
is the initial extended tensor product mesh). As shown in the proof of Lemma \A.2\ 
if there is no isolated subdomain at level i — 1, then one can order the l-edges in E* 
such that 

^ ^ f m — 1, if ii^j is a horizontal l-edge 
^ ~ \^ n — 1, if ii^j is a vertical l-edge 

If a subdomain ^ is isolated, then one can order the l-edges in the subdomain ^ 
such that the above equality holds for all but one vertical l-edge in S^. For the 
exceptional l-edge, 

Vj' = n — 2. 

Figure [9], illustrates an example for the order of l-edges in E^ in the case of 
m = 4, 77, = 3. 

From the above facts, one has 
E dimW[Ej] = E - ^ + 1) + E -n+l) + 5.,, i > 1, (18) 

B.GE' EjGEi E,€El 

where E\ is the set of horizontal l-edges in E*, E* is the set of vertical l-edges in E*, 
and 6i is the number of isolated subdomains at level i — For i = 0, since E*^ is the 
set of l-edges of a tensor product mesh it is direct to check that 

E dim W[Ej] = E (^Z - ^ + 1) + E - + 1) + (m - l)(n - 1). (19) 

EjGEO EjGEl EjGEO 

Now the dimension formula dUD follows from ([I3]), (II7D, ([T8]) and ([19]). □ 
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4.3 Examples 

In this subsection, we give two examples to count the dimensions of the sphne spaces 
S(m, n,m - l,n - 1, 5^ „). 

Example 4.1. In the case o/m = 2,n = 2, ^m,n is a regular hierarchical T-mesh 
and we denote it as ^ . By the dimension formula f[T^ . 

dim S(2, 2, 1, 1, ^) = dim S(2, 2, 1, 1, 

= V+ - E + 6 + 1 

where, is the number of crossing vertices, E is the number of interior I- edges, 
and 6 is the number of isolated cells of S/^'^ . This formula has been presented in j^. 

Next we give a concrete example for the case m = 3,n = 3. 

Example 4.2. In this example, the hierarchical T-mesh S^-^^j, (which is the same as 
the T-mesh in Figure^c)) and its extension are illustrated in Figure\^ There 



^3,3 ^3,3 

Figure 8: A hierarchical T-mesh ^ 3 and its extension. 

are three isolated subdomains (one at level and two at level 1) in the middle of 
Thus 5 = 3. It is easy to count that = 166, Eh = 21 and Ey = 19. So 

dim 8(3,2,2,2,5^3^3) 

= V^-{3- 1)Eh - (3 - l)Ey + (5 + (3 - 1)(3 - 1) 
= 93. 

We can also get the dimension by counting the number of basis functions. 

There are (12 — 4) (13 — 4) = 72 basis functions over the 11 x 12 tensor product 
mesh =^^3. At level 1, 7 basis functions are added to the basis functions set by 
subdividing two subdomains. At level 2, 14 basis functions are added to the basis 
functions set by subdividing three subdomains. Totally, there are 72 + 7 + 14 = 93 
basis functions, and therefore the dimension of S(3, 2, 2, 2, ^■^^ is 93. 
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5 Conclusions and future work 



This paper presents a dimension formula for the sphne space S{m,n,m — l,n — 
1, over certain type of hierarchical T- meshes. By using the smoothing co-factor 
method, we set up a bijection between the spline space S(m,n, m — l,n — 1, ^) 
and a conformality vector space W [£/']. Then we decompose W[^] into direct sum 
of simple linear spaces by using homological technique, and a dimension formula is 
thus obtained. At a by-product, we also obtain a set of basis functions for the spline 
space S(m, n, m — 1, n — 1, ^7). 

It is important to construct a set of proper basis functions which have nice prop- 
erties for the spline space in applications. We will study this topic in the future. 

6 Appendix 

In this section, the proof of Lemma 14.21 is given. 

Proof. We first outline the proof strategy. For getting the result, we choose a set of 
bases {k*} for and construct a function Bj{x,y) G S(m, n, m — 1, n — 1, ^)(= 

W[^]) such that the conformality vector of Bj{x,y) corresponding to the vertices 
of Ei is k;.. Then vr : W[£^] — > W[Ei] is surjective. 

Before we discuss the construction of Bj{x,y), some preparations are given as 
follows. 

Preparations 

1. Choose a set of bases for iy[£'j]. If Ei is horizontal, let A^. : — oo < a^i < 
X2 < ■ ■ ■ < Xr^ < OO be the x-coordinates of the vertices of Ei on ^i. The 
B-spline basis functions {N[xj,Xj+i, - ■ ■ , Xj+m+iW/Ji^''^ is a set of bases of 
S(m, m — l,Ax) and {kjYj'Sj^'^ is a set of bases of where k* is the 
conformality vector of N[xj,Xj+i, ■ ■ ■ ,Xj+rn+i]- For a vertical interior 1-edge, 
we can choose a set of bases of similarly. 

2. Define an operator on T-meshes. The operator D on a T-mesh ^ associated 
with an l-edge Ei is defined as follows. D£;-(e^F) = (Ei, where is the T- 
mesh by removing E^ from Then ttioD^.^^/') = Ei , and 7r2oD£;-(^) = 

3. Define an order of E. The interior l-edges E of ^ „ can be divided into disjoint 
sets E\i = 0,l,...,l, i.e., E = U[^qE' with E^ n E^' = 0, z ^ j. is the set 
of interior l-edges of ^ „ defined at level i, i = 0,1, . . . ,1, where / is the level 
of Xi,n- The order 'between E^ and E"? is defined as: Ve^ G E^, Ve« G E", 
a q < p, then e'^ -< e^. In the following, the order "-<" inside E^ is presented. 
Consider E^, and suppose the T-mesh is obtained by deleting the l-edges 
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U-^p^^E' from Without loss of generality, it is enough to consider the 

case of p = Z. 

Then, the order inside E' is defined as follows with respect to the case of Z = 
and / > 0, respectively. 

• / = 0. The T-mesh =^ is a tensor product mesh The progress of 
choosing interior l-edges sequence is presented in the following. 

Se ^ the set of interior l-edges of ^ 
i ^ 1 

while 7^ do 

if there is no trivial interior l-edge then 

choose any interior l-edge 
else 

choose a trivial interior l-edge (i.e Vr[-E'?] = 0) 
end if 

^ ^ 712 0-0 Efi^) 

Se <— the set of interior l-edges of ^ 
i i + 1 
end while 

By this progress, we get a sequence of E° denoted as E^, E2, ■ ■ ■ , where 
E^ is chosen earlier than E^^^. The order is defined as -< E^ 

• I > 0. Select an l-edge £^ e E'. Then E must intersect with an (m,n)- 
subdomain D at level (Z — 1). If i? is vertical, we consider all the vertical 
l-edges in E' which intersect with D and arrange them from left to right. 
Then E is the i-th l-edge. By the structure of E is also the i-th 
l-edge in another (m, n)-subdomain at level (/ — 1) which intersects with 
E. Hence we can denote the position of as L{E) — i. Considering the 
local tensor product structure on D at level — 1), L(E) — i < m. If £" is 
horizontal, L(E) can be defined in a similar way by arranging horizontal 
l-edges in E' from top to bottom and L{E) < n. 

Define l-edge sets Ai, A2, A3, A4, A^ as 

Ai ^{E G E' : If ^ is vertical, L{E) < m - 2; 

If E is horizontal, L{E) < n-2}, 
A2 ^{E e E' : £; is horizontal, L{E) = n - 2}. 
A3 ^{E e E' : £; is vertical, L{E) = m - 2}, 
Ai ^{E e E' : £; is horizontal, L{E) = n - 1}, 
A^ ^{E e E' : £; is vertical, L{E) = m - 1}. 

The progress of choosing interior l-edge sequence is presented in the fol- 
lowing. 



17 



forJ^= 1 to 5 do 

^ the set of interior 1-edges of ^ in 

while 7^ do 

if there is no trivial 1-edge in Se then 

choose any interior 1-edge E^j from Se 
else 

choose any trivial 1-edge Elj from Se 
end if ^ 

^^VTaoD^. (^) 

S^; ^ the set of interior l-edges of ^ in Ai 

J ^ J + 1 
end while 

end for 

This progress generates an interior 1-edge sequence E[ -^^, E{ 2,- ■ ■ , E[ j^, 
-^2,15 -^2,2) -^2,j2' ^5,35- 1-edge ei is chosen earlier than 

another 1-edge 62 in the sequence, we define 62 -< ei. 

For an example, see Figure [9l we consider the case of / = 2. Here m = 
4, n = 3,p = 6, g = 6. An 1-edge in the T-mesh belongs to the set in the 
right side of Figure M whose name has the same color as the 1-edge. And 
its order of an interior 1-edge sequence of is labeled near itself. 
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Figure 9: An example of E^, where m = 4, n = 3, p = 6, g = 6. 



BUx,y) construction 
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Now we are ready to construct Bj{x,y). 

When / > 0, the B-sphne function Nj{x,y) can be constructed by the tensor 
product mesh over (m, n)-subdomains at level (/ — 1) contained e. Here Nj{x,y) is 
related to the knots of an interior l-edge e in the sequence associated with y4i(or 
^2,^3). It is a tensor-product B-spline function, which equals to Bj{x,y) up to a 
nonzero constant. For trivial interior l-edges, no function is introduced. Suppose 
that we construct the functions associated with all the l-edges whose orders are 
"larger" than the order of E in the sequence associated with A4. We now construct 
the functions associated with E. E is a non-trivial horizontal l-edges and we denote 
the x-coordinates of E as {xi, X2, - ■ ■ , Xd}- Let a be the number of vertices formed 
by two l-edges with level / intersection. This type of vertex is called an (/, /)- vertex 
later. The i-th knots sequence of E's x-coordinates is {xi, ■ ■ ■ , Xi+^m+i)}- According 
to the structure of the current mesh, a G {0, 1, 2}. 

1. a = 0. If Z > 2, a B-spline can be constructed over a tensor product mesh over 
all the (m, n)-subdomains at level (/ — 2) containing E^s i-th knot sequence. 
If / = 1, we can construct a B-spline according to the tensor product mesh of 
level 0. The required function Bj{x, y) is the same as the B-spline function up 
to a nonzero constant. 

2. a = 1. If Z > 2, consider the unique (/,/)— vertex, denoted as Q. Assume 
the (m, n)-subdomain which contains Q at level (/ — 2) is S. If the interior 
vertical l-edge that contains Q goes through S, we can construct a B-spline 
similar to the case of a = 0. If it does not go through E, one of its endpoint 
must be within S and another is outside of S according to the structure of 
the special hierarchical T-mesh. Suppose the intersection point between the 
vertical l-edge and S is R. We extend the l-edge from the endpoint within 
E and cut the boundary of S at P. Then the n + 2 vertices can be used 
to construct a B-spline Ni{x,y) by a tensor product mesh over all {m,n)- 
subdomains containing the given (m + 2) vertices. These n + 2 vertices consist 
of P, R, and n vertices in the original vertical l-edge within S. Therefore, 
the conformality vector of Ni{x,y) is nontrivial at these (m + 2) vertices and 
trivial at the new vertices without P in the process of extending. Suppose 
ki is the conformality factor at P. If another function g{x,y) is constructed 
such that its conformality factor at P is —ki and the conformality factors at 
the new vertices without P and the given (m + 2) vertices are zeros, then 
f{x,y) = Ni{x,y) + g{x,y) differs by a multiplier A (7^ 0) with the function 
Bj {x, y) which we need by the linear property of Jt^ and Example 13.21 

Now we discuss the construction of g{x,y). If we remove Q from the set 
consisting of P and all the vertices in the original vertical l-edge, the number 
of remaining vertices is not less than (n + 2). We choose 2)— vertices from 
these remaining vertices and P must be one of them. Similar to the case of 
a = 0, there is no (/, /)-vertex in these (n + 2)— vertices. So, a B-spline N2{x, y) 
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is constructed and we denote ^2 as the conformality factor of N2{x,y) at P. 
By Example 13.21 k2 ^ and the conformahty vector of N2{x,y) is trivial at 
the new vertices without P and those given (m + 2) vertices according to Q 
removed. Then we define g{x,y) = —ki/k2N2{x,y). 

For an example, see Figure [TOl where m = 3, n = 3. Here P is shown as "□" 
in the figure. Ni{x,y) is constructed over the tensor product mesh shown in 
the upper right of this figure and the one in the lower right is used to construct 
N2{x, y). Compute ki of A''i(a;, y) and k2 of N2{x, y) at P as we have shown in 
Examples 13.11 and 13.21 Then f{x,y) is given by 

f{x,y) = N,{x,y)-'^N2{x,y). (20) 
k2 

When / = 1, a B-spline associated with the given m + 2 vertices can be 




Figure 10: Construction of the linear combination of B-splines in the case of m 
n = 3. 



constructed by the tensor product mesh of level 0. 

3. a = 2. If / > 2, there are two (/, /)— vertices Pi and P2 in the given (m + 2) 
vertices. One of the two vertical l-edges at Pi, P2 respectively must go through 
S' that occupied by (m, n)-subdomains containing Pi and P2 at level (/ — 2). 
So, we can construct the function by the method presented in the case of 
a = 1. When / = 1, the function can be constructed by the tensor product 
mesh of level and the case I > 2. 
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When all the 1-edges in the sequence associated with A4 are deleted, there arc no 
(/, /)— vertices at the 1-edges in A^. So we can construct the function we need by the 
method presented in the case of a = 0. 

When / = 0, by the order defined on E° and the tensor product mesh at level 0, 
B-spline surfaces can always be constructed. So as stated previously, tt is surjective. 

□ 
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